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Abstract — We prove that a group of injection-locked oscillators, each 
modelled using a nonlinear phase macromodel, responds as a single 
oscillator to small external perturbations. More precisely, we show-that 
any grDun-pf injection-locked oscillators has a single effective PRC ||l|] or 
PPV [H], [Hp that characterises its phase/timing response to small external 
perturbations. This result constitutes a foundation for understanding 
and predicting synchronization/timing hierarchically in large, complex 
systems that arise in nature and engineering. 



I. PRC/PPV Phase Macromodels 
Given an ODE or DAE description 



dt 



q(x(t)) +m+b(t)=Q 



(1) 



of an oscillator with an orbitallv stable T-periodic autonomous 
solution x s (t), it can be shown jg], [^] that the timing jitter or 
phase characteristics of the oscillator, under the influence of small 
perturbations b(t), can be captured by the nonlinear scalar differential 
equation 



dt 



a(t) = vj(t + a(t))-b(t), 



(2) 



where the quantity vi( ), a r -periodic function of time, is known 
as the Phase Response Curve (PRC) [|l]] or Perturbation Projection 
Vector (PPV) 

For convenience, we scale the time axis to normalize all periods to 
1. Define a 1 -periodic version of the steady state solution to be 



Xp{t)=x s {tT), 
and a 1 -periodic version of the PPV to be 
p(t) = v 1 (tT). 



(3) 



(4) 



Using these 1 -periodic quantities and defining f — j, (Q) can be 
expressed as 



dt 

Defining phase to be 



o(0 = p T (ft + fa(t))-b(t). 



<H0 = ft+fa(t), 



(5) 



(6) 



(^) becomes 



j t *<!) =f+fp T (Ht))-b(t)- 



x(t), the solution of (|l]), can often be approximated usefully by a 
phase-shifted version of its unperturbed periodic solution, i.e., 



<^(7>(O)=%(0(O). 



(8) 



(g) (or equivalently, (|Vp) is termed the PPV equation or PPV phase 
macromodel. In the absence of any perturbation b(t), note that a(t) 
(w.l.o.g), x(t) = x s (t)=x p (ft) and <j)(t)=ft.We will call the latter 
the phase of natural oscillation and denote it by — fi- 
ll. Derivation of hierarchical PPV macromodel 
A. Coupled Phase System and its Properties 

1) Coupled system of PPV phase macromodels: Consider a group 
of N > 2 coupled oscillators (Figure |3]). We model each oscillator by 
its PPV equation (0): 



dt 



fc(0 =fi+fi?i(<i>i(t))-b i (t), / = !,• 



(9) 
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Fig. 1. Oscillator system with internal coupling and external inputs. 



where i-subscripted quantities refer to the i th oscillator. Inputs to each 
oscillator are drawn from two sources (as depicted in Figure [l^: 1) 
internal couplings with other oscillators, and 2) external sources. bi(t) 
can therefore be written as 



-a, t 



+ I h («>y(0) , 

7=1 



(10) 



where 3,(0 is the external input (i.e., from outside the group of N 
oscillators) to the ! th oscillator, and b/j (<j>j(t)) represents the influence 
of the / h oscillator on the ! th . 

We make the natural assumption that the bu(-) are 1 -periodic — 
i.e., that each oscillator generates outputs that follow its own phase 
and timing properties; it is these outputs that couple internally to the 
inputs of other oscillators. Note that as i varies, the dimensions of 
Pi(t), 3j(t) and bjj can differ, since they depend on the size of the 
( th oscillator's differential equations. 

The system of N equations (^) can be written in vector ODE form 
as 

(11) 



dt 



0(0 =2*(0(O)+M#(O 



(7) where 



?(0 = 



and b*($,t) = 



>l(0 
>jv(0. 

" /l+/iPf(0i) -L^iblj^M 

fN+f N p T N (M ■ L^ 1 b Nj (*j(0) 
7ip[(0i)-ai (0 

fNP T N (<l>N)-a N (t)_ 



(12) 



(13) 



(14) 



We will refer to ( |1 1|) as the Coupled Phase System ( CPS). 
Note that (0,0 = in the absence of inputs external to the group of 
oscillators, i.e., when a/(t) =0. Note also that gd(<j>) is 1-periodic in 
each component of 0, i.e., it is 1-periodic in each 0,-. Such functions 



2 



are termed cylindrical ^ page 236]. 

2 ) Locking in the absence of external inputs: 

Assumption 2.1: In the absence of external inputs (i.e., if 5,-(f) = 
0,V0, assume that the group of N oscillators locks stably^ to a 
common frequency /* (equivalently, with a common period T* = A). 
Denote the phase of the ; th oscillator, when locked as in Assumption 
Ej] to the other oscillators in the group, by <j)*(t). Note that this 
phase will typically be different from the oscillator's phase of natural 
oscillation <j)f(t) = ft, on account of inputs via coupling from other 
oscillators in the group. 
Denoting 

W(01 



0*(O = 



KO. 



note that 



"(f) satisfies ( |ll[ ) with no external inputs, i.e., 
d 



dt 



(15) 



(16) 



We term <p*(t) the system phase during externally-unperturbed lock. 
3) D-periodicity of <j>*(t): T*-periodicity of each oscillator locked 
at frequency /* implies that 

is r*-periodic Vi; i.e., 

Xi(t + T*)=xi(t). or 

x P M('+T*))=x L P .M(t))- <IS ' 

From definition, x p j(-) is 1-periodic. ( fl8| ) is satisfied for arbitrary 
1-periodic x p j(-) iff 

^*{t + T*)=n + ^(t), Vf, neZ. (19) 

Define the ideal phase of oscillation at frequency f * to be 

&*l{t)=ft. (20) 

0ideal(O sa ti snes ( p^l ) (with n — 1), but note that the phase of each 
locked oscillator in the system need not necessarily equal 0;d ea l(O- 
A more general form for (j>*(t) that also satisfies (|19j) is 

0*(O = feW+A^W = f*t+Atf(t), (21) 
where A0*(?) is itself 7"* -periodic^} Equivalentljj^J 

" (O = 0idealW+A?*(O=.n+A?*W, (22) 

where 

rAtf(*)l 



A0*(O = 



(23) 



A0*(r)_ 

is r*-periodic. Functions of the form ( ^2[ ) are termed D-periodic or 
derivo-periodic with period T* [||]. {A0*(j)} represent short-term 
phase changes within each cycle that do not affect the long-term 
frequency of the oscillator. 
The above considerations motivate: 

Assumption 2.2: $*(/), the phase of the CPS during externally- 
unperturbed lock, is r*-D-periodic. 

4) Arbitrary time shifts of <j)*(t) are also solutions of the CPS: 
Lemma 2.1: The phase during externally-unperturbed lock, <j)*(t), is 
not unique. Indeed, for any arbitrary time-shift T, 



(24) 



1 Section G-E will expand on the notion of lock stability. 
2 A0*(f) can in fact itself satisfy ( |19| ) with arbitrary n, but n j£ would make 
the long-term frequency of Xp,i($j(tY) different from /*, violating Assumption 



We use the notation, borrowed from MATLAB, that the sum of a scalar 
and a vector means that the scalar is added to each element of the vector. 



solves (|1J). 

Proof: Follows directly from substituting ( [m| ) in the autonomous 
system ( |l6| ) and using the facts that 1) g^(-) in ( |l3| ) is cylindrical 
with period 1, and 2) <j)*(t) is r*-D-periodic (Assumption 2.2|. ■ 



B. Periodic time-varying linearization of the CPS 

1) Linearization under small-deviation assumption: If the external 
inputs {ai(t)} are small, then t) is small and (JuJ) constitutes a 
small perturbation of (|l6|). We express <j)(t), the solution of (|ll|), as 
a deviation from <j>*{t), the solution of (flfl): 



(25) 



We term 8tj)(t) the orbital deviation. Using (|25j), we now attempt to 
solve ( p~jj ) via linearization. 

Assumption 2.3: 8tj)(t) remains small for all t, provided the external 
input Sa(v) is small enough for all r.n 



Applying Assumption 2.3 , we start the process of linearizing (|1 1|): 

d 



iH y(t)+-sHt) ^gt{r(t))+^(r(t))m) 



(26) 



Using (Jlfj), we obtain 



-«0(O ^-#(**(O)«0(O+**(*(O>O 
*^(?*(t))**(0+^(**(0,0 

dip 



= (^(?*W) + ^(?*(0,0)^(0 



(27) 



+ 



■(0,0 



From the definition of 

%{■,•) 0. observe that ^|(0*(f),O is a 
diagonal matrix with entries 

i.e., it is directly proportional to the external inputs {«,■(?)}, which are 

small by assumption. Therefore, the product term ^i- (0*(f),f) 5<j)(t) 

in ( ^7| ) is of second order and can be dropped from the linearization. 
Applying this observation and denoting 



( ^7| ) becomes 



/;(O = 4f(0*(O,O, and 
Sext(0 = fy(?*(0>0> 



50(f) ~^(f)50(O+iext(O 



(28) 
(29) 

(30) 



( po| ) is the linearization of the CPS ( [TT] ) around its externally- 
unperturbed solution <j)*(t). 

2) T* -periodicity of J^(t): From (|T^), we can obtain expressions for 
the entries of 7| (f). The diagonal entries of Jl are 



j; u (t)=fii? t T w(t)). j>,- (*;«), 

4 i.e., < M||2ty(-, -)|| for some finite constant M > 0. 



(31) 
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while the off-diagonal entries are 



(32) 



Because of the 1 -periodicity of /?,(■) and bjj(-), and the T*-D- 
periodicity of (j>*(t), each entry of is r*-periodic, hence the entire 
matrix function J^(t) is T* -periodic. The linearized CPS ( |3C| ) is 
therefore periodically time varying, i.e., it is a linear periodically 
time varying (LPTV) system. 



can be rewritten using v,(-) and «,•(■) as 

AT 

**(0 = L%(0^ (f ""' b) iffo) •**<>+ 



scalar 



(43) 



V- (T)-fcext(T) dT. 



scalar 



C. T* -periodic homogeneous solution of the linearized CPS 
Lemma 2.2: The homogeneous part of the linearized CPS ([?(]), i.e., 

j t H{t) =4(t)S$(t), (33) 
has the T* -periodic solution 



(34) 



Proof: Follows immediately from differentiating (|l6| ). Note that 
7"*-D-periodicity of <j>*(t) immediately implies that 8<p*(t) in (Q) 
is r*-periodic, since 



8**(f)=/* + ^A**(0, 



(35) 



with the latter term r*-periodic. 



£>. Floquet-theoretic solution of the linearized CPS 

Floquet theory [Q] provides an analytical forrrj^ for the solution of 
(§): 

5^(t) = U(t)D(t-t )V T (t )5lo+ 

f t -* (36) 

U(t) / D(t-T)V T (t)h ext (T)dT. 
J to 

U(t) and V T (t) are r*-periodic matrix functions, of size NxN, that 
satisfy 

(37) 



U(t)V T (t) = V T (t)U(t)=I NxN . 



( p7| ) implies that the columns of U and V are bi-orthogonal, i.e., 



vJ(t)-Uj(t) = Sij, i,j=l, 



,N. 



( p7| ) can be written more explicitly, showing V; and 5j, as 

/ : : ^ /I 

«i(0 ••• m(t) 



\ 



j 



\ 



Note that, in particular, 

v[(f)-Si(r) = l,Vf, 

a relation we will rely on later. 

£>(■) in (^) is a diagonal matrix of the form 

D(t) = 



(38) 



(39) 



(40) 



(41) 



V 



where {p,} are termed Floquet or characteristic exponents, and 

p,.A e Ai,r 5 ,-=1,...,^ (42) 

are known as Floquet or characteristic multipliers. Note that D(f) is 
not periodic. 

5 ( pij ) holds for the case where the Floquet multipliers (jllj) are distinct, but 
subject to Assumption E.4, all subsequent results in this section remain valid 
for the case of repeated Floquet multipliers. 



E. Conditions on Floquet multipliers; stability and isolation of<j>*(t) 



Lemma 2.1 



suggests that the stability of the CPS ( |11| ) around its 
locked solution (|15j) in the absence of external inputs is of an orbital 
nature [||, Definition 5.1.1] and not, e.g., of a Lyapunov one |^| 
Definition 1.4.1]. 

Lemma 2.3: At least one of the Floquet multipliers {p,} ( p2[ ) is 1 
(equivalently, at least one of the Floquet exponent {/i,} is 0). 

Proof: Follows from the existence of the T*-periodic homoge- 
neous solution of the linearized CPS ([35]), established in Lemma 2.2 , 
and [§ Corollary 2.2.3]. ■ 
Lemma 2.4: |p,| < l,Vi (equivalently, M.(u.f ) < 0). 



2.1 



(i.e., that the externally- 



Proof: Follows from Assumption 
unperturbed oscillator system is mutually injection locked in a stable 
orbit) and |^| Theorem 5.1.3], ■ 
We now make an additional technical assumption regarding the 
Floquet multipliers: 

Assumption 2.4: Only one Floquet multiplier (pi, without loss of 
generality) is 1 (equivalently, w.l.o.g, pi = 0). 



There are several 



Note that Assumption strengthens Lemma 2.3 . 

factors that motivate this assumption: 

1) Assumption 2.4, together with Lemma 2.4, constitute sufficient 
conditions for the Andronov-Witt theorem [^, Theorems 5.3.8 
and 5.1.2] to hold. The Andronov-Witt theorem guarantees 
that the CPS ( p"i| ) is asymptotically orbitally stable (a.o.s) [5, 
Definition 5.1.2] with the asymptotic phase property (a.o.p) [5, 
Definition 5.1.3]. These properties are central to the intuitive 
concept of stable lock, assumed in Assumption 2.1 and typi- 
cally valid in most applications. 

Assumption |2.4| also constitutes a sufficient condition for the 
orbit (j>*(t) (|15|) to be isolated, i.e., not embedded in a contin- 
uum of orbits with continuously-varying periods Theorems 
5.3.9 and 5.2.3]. Isolation is also central to the intuitive notion 
of stable lock. 



2) 



3) 



Although Assumption 2.4 is not a necessary condition for 
asymptotic orbital stability with the asymptotic phase property, 
or for isolation, oscillators that do not satisfy the assumption 
while still being a.o.s + a.o.p tend to be "unnatural". For 
example, |^| Example 5.2.2] features an orbit that is a.o.s + 
a.o.p but has three repeated Floquet multipliers that equal 1; 
however, this orbit is not periodic. 



Thus, in most practical situations, Assumption 2.4 is essentially 



equivalent to the oscillator group's being "stably locked". 



F. Assumption 2.3 (deviations are small) is invalid 
Lemma 2.5: Assumption 2.3 is invalid. 



Proof: Using Assumption 2.4, the second summation term of 
which captures the linearized system's response to external 
perturbations b ex t(t), can be written as 



scalar c\ (t) 



Ui(t) [ V 
J t 



+ ZUt) t e^>-^vJ(T)-b ext (T)dT. 

i=2 



(44) 



scalar c;(f) 
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The first term is of the form ci(f)f?i(f), where c\(f) is the scalar 

rt 



to 



(-r)-feext(-r) dr. 



Because vi(t) is periodic, there exist many possibilities for small 
^ext(0 that make ci(f) increase indefinitely and without bound as 
f increases. For example, if f5 ext (f) = £u\(t), with e ^ being 
any constant, then, from (|40|), ei(f) = (f — jn)e; ci(f) increases 
without bound. The remaining terms in ( |44| ) are bounded because 
R(/x,) < 0, Vf > 1, hence cannot cancel the unbounded increase of 
the first term. 



In other words, Assumption 2.3, upon which the linearized system 
(f3(j|), its solution (|3rf), and indeed, the expression for the unbounded 
term c\(t) in ( p3| ) all depend, is violated. Thus we have arrived at 
a contradiction, implying that the original premise Assumption 2.3 
must be invalid (subject to the other assumptions' validity). ■ 
That deviations can grow to be large even when external inputs 
remain small is a manifestation of the inherently marginal nature 
of orbital stability, i.e., that /ij =0. 



G. Time-shifted perturbed response assumption 
Lemma 2.6: Without loss of generality, 



Proof: The first summation term in (|43|), i.e., 



(46) 



(47) 



represents a general solution of (|33|). We already know that 50* (f) 
( p4| ) is a nontrivial periodic solution of (^). Using p\ = 1 from 
Assumption 2.4, this periodic solution must equal the i = 1 term 
in (pj|), since (also from Assumption |2.4j ) the remaining terms for 
(' = 2, ■ • • ,N are not periodic and indeed, decay to as t — > °°. Hence 
we have 



{t)=k z u\ (t)v\(t )-8(f>*{to). 



(48) 



scalar constant k\ 



where ki is an arbitrary scalar constant. Note that k\ ^ 0, otherwise 

§ would be identically zero, hence would not match any nontrivial 
(?). Choosing ki = i- (without loss of generality, since v\ (t) can 
be scaled to satisfy (^)) results in Q). ■ 
Geometrically, 4i^*(t) is the tangent to the externally-unperturbed 
orbit of the locked system in phase space; Lemma 2.6 thus justifies 
the terminology tangent vector for u\{i). 

Attempting to restore validity to the failed linearization procedure 
above, observe that using G46j), the unbounded term in (|45|) can be 
written as 



: C1 (o-rw- 



ci(0«i(0 

Observe also that if c\ (t) were bounded and small, then 

$*(t)+ci{t)-r$*(t) ~$*(t + ci(t)), 



(49) 



(50) 



to first order. This suggests that the unboundedly growing component 
of 8(j)(t) in @ may be the manifestation of a time-shift to the 
unperturbed solution 0*(f). A time shift alone the orbit is also 
suggested by the definition of orbital stability pi Definition 5.1.1] 
and by the physical intuition that autonomous oscillators, having no 
intrinsic "time reference", can "slip in phase", i.e., they cannot correct 
errors in phase. Accordingly, we modify the assumed form of the 
perturbed solution ( |25| ) to 
Assumption 2.5: 



where 8<j)(t) remains small for all time (i.e., |50(f)|| < M||fo e xt(f)|| 
for some finite M > 0). 



(45) 

cc g (t) is a (yet-to-be-determined) time shift that can depend on the 
input perturbation b ext (t) and can grow u nbo undedly with time. 
Importantly, we have retained Assumption 2.3 i.e., that S(j)(t) in 
Slh remains bounded and small for all time. 



We shall prove that unlike ([25]), the time-shifted deviation form ( |5l| ) 
will allow 50(f) to remain bounded and small, providing the time- 
shift a g (t) is chosen appropriately. 

H. Base for time-shifted linearization 



In Section II-B, the CPS was linearized around the unperturbed 



orbit 0*(f). The process of linearization relied on the fact that 0*(f) 
satisfied (|l6[). We would like to find a replacement for (JT^) that is 
satisfied by 

0(f) = 0*(f + a g (f)) (52) 

instead. 

Lemma 2.7: Given any scalar, differentiable function oc g (t), the CPS 
( p"l| ) is solved exactly by ( p^ ) for perturbations of the form 

bt($(t),t)=K(t)r n (t + a g (t)), (53) 

where K{t) = f t a g {t). 

Proof: Denoting "shifted time" to be 

ft 4 * + (%(»), (54) 

substituting ( |5^ ) and (p3|) into (|ll|) and simplifying using ( |l^ ) and 
(fill, we obtain 

d r» 



(l + a g (t))^ 0*(ft) =^(0*(ft))+X(f)w 1 (ft; 
K(t)»i(S) 



■ ^ d 
>a ^diT 



(55) 



=^%(f)«i(ft) = iT(f)«i(ft) 
^>a g (t)u] (t + a g (t)) = K{t)u\ (t + <x g (t)). 
55]) is always satisfied if cc g (t) and K(t) are related by 
h g {t)=K[t). 



(56) 



/. Time-shifted linearization 

We proceed to linearize the CPS ( pi] ) around solutions of the form 
To this end, we split the external input S^(0,f) ( |l4"| ) into two 

parts: 

b^,t)=b H {t) + b^ 2 (4,t), (57) 

with the intent that if only the first component b§ . (t) is retained, 
then @ should solve the CPS ( Jl~T| ) exactly, i.e., 



j t <t>*{t + a g (t)) 



i*(t+a g (t)))+b Pl {t). 



(58) 



(51) 



Motivated by Lemma 2.7, we explore perturbations of the form ( |53[ ) 
along the tangent vector, i.e., of the form 

b^(t) = K(t)ui(t + a g {t)). (59) 

Given any small external perturbation b ^ (0 (f), f) (|l4|), our goal is to 
find such an a g (t) (and, using Lemma |2.7[ its derivative K(t)) that 
S^ 2 (',') m (plj), as well as the orbital deviation 50(f) in (p"i|), both 
remain small. 

The flow of the time-shifted linearization procedure is: 

1) Start by assuming any scalar function A*(f); 

2) Define a g {t) using ( p6|), i.e., £ a g (t) =K(t); 

3) Define o^j(f) using (|59|); 

4) Incorporate the split-up form (^) of the external input pertur- 
bation in the CPS (0); 
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5) Assuming a solution of the form ( j5l| ) in Assumption 2.5 
linearize (|llj) using ( pi) ) as the base case; and 

6) Using the solution of the above linearization, obtain a constraint 
on a g (t) (equivalently, on K(t)) under which Assumption 



2.5 



holds with 50(0 bounded and small for all time. The 
equation specifying this constraint will turn out to be of central 
importance, in that it governs the phase/timing responses of the 
injection-locked system of oscillators to external perturbations. 
7) When oc g (t) (equivalently, its derivative K(t)) is chosen to 
satisfy the above constraint, show that the phase deviation 
MO in Q always remains small, thus validating Assumption 



2.5 



and the entire time-shifted linearization procedure. 
Starting from Step Q write the CPS (JTT|> as 

jj{t) =2f(?(0)+^j(0+Zfe(?(0. 

Incorporating (|l]), @ and Q in we obtain 



d rr 

dt 



$*(t + a g (t)) + S<l>(t)] =g <) ($*(t + a g (t)) + 8<l>(t)) 



(60) 



(61) 



+K(t)u 1 (t + a g {t))+b$ 2 + a g (t)) + 50(0,*) 
Linearizing g^(-) in ([5l|), we obtain 

- $*(t + a g (t))+- 50(0 = ft + «%(*))) 
+ J$(t + a g (t)) 8<j>(t)+K(t)ui (t + a g (t)) 

Applying the base for time-shifted linearization ( p8[ ) and our proposed 
form ( |5^ ) for b^^t), ( jo2| ) can be simplified to 

d 



(62) 



dt 



m) 



-a g (t))d^(t)+b^(r(t + a g (t)) + 8<l)(t),t). (63) 



Observe that from definition (|57|), (p9|), 

fy 2 (0* (t + a g (t) ) + 50 (0 , t) =b$ (0* (f + a g (t) ) + 50 (t ) , t) 



-K{t)ui(t + a g {t)), 



(64) 



hence ( |63| ) can be written as 



50(0 



(65) 



a g {t)) 50 (0 - K{t)u\ (t + a g (t)) 
Using the same reasoning as for (§]) in Section |lI-Bl| , 50(0 

in the 

last term of (^) can be dropped because it contributes only a second- 
order term to the linearization. Hence (p5D becomes 



- 50(0 =4(t + <%(*)) **(0 -K{t)m (t + <%(*)) 



+ fy,(0*(r + < % (O),O 



(66) 



= 7;( ? r )50(o+^(0*(f l "),O-^(O«i(' 1 ") 
= y;( f t )50(o+^ 2 (0*(/ t ),O, 

where we have used the notation f , defined in (154), for shifted time 



J. Recasting time-shifted linearization in LPTV form 

We would now like to obtain an analytical solution of (^) and use 
it to validate that 50(0 remains small for all time. However, two 
differences bet ween ( pOp and (^(j) make this more involved than for 
@ in Section [il-Pf 

1) The input to ( |66[ ) is b^ 2 (-,-), not b^,(-,-) as in (pp[). Whereas 
the latter is known small (due to the assumption of small 
external perturbations {a;(0} m dl^X there is no guarantee 
that £>0 2 (-,-) is also small. Ensuring that b^A-,-) is small is 
important: if even the input to (p6|) cannot be guaranteed small, 



it is unreasonable to expect that its solution will remain small 
for all time. 

2) Unlike @, which is LPTV, (|3) is not LPTV because though 
7^(0 is r*-periodic, J*^(t + ag(t)) is not, except for special 
choices such as a g (t) = 0. We are interested in a solution of 
(p6|) that is valid for any <X g (t) (equivalently, any K(t)), if 



possible. Becau se 
in Section 



is not LPTV, the Floquet expressions 



1I-D do not apply directly. 



Both issues can be addressed by restricting a f (0 = K(t) to be small. 
We state this as an assumption for the momentrt 
Assumption 2.6: 

K(t)±a g (t) (67) 



is small and bounded with respect to b§( 
particular, \K(t) \ < 1. 



(Q for all time. In 



The first consequence of Assumption 2.6 is that it becomes possible 
to guarantee that £a>(v)i me m P ut t0 me time-shifted linearization 
(pq), is small: 

Lemma 2.8: S ( (, 2 (0*(? t ),O is small for all t\t. 
Proof: From definition (p7[), 

b H {4*(?),t) = ^(0*(/ f ),O -K{t)uitf). (68) 

The first term is small from our underlying assumption of small 
external perturbations. The tangent vector ti\ (t) is a periodic, bounded 
quantity, hence under Assumption 2.6, the second term is also small. 



Another important consequence of Assumption 2.6 is 
Lemma 2.9: The mapping ( ^4| ) 

/i-^r"'", i.e., f f (0 =t + a g (t) 

is invertible. 

Proof: It suffices to show that the mapping is a monotonically 
increasing one, i.e., its derivative is always positive. We have 

d 



dt 



■ft 



1 



■a g (t) = l+K(t). 



From Assumption 2.6, \K(t)\ < 1, hence 4;fi(t) > 0, i.e.. 

monotonically increasing. 

We now make the following definitions: 

50V) ^50(0, 



b\ 2 (a,t T ) =b^ 2 (a,t), and 
b\{aj) =b${a,t). 



v(t) is 
■ 

(69) 
(70) 
(71) 



The significance of the invertibility of shifted time t 1 , as established 
by Lemma |2.9j , lies in that the above definitions become possible: 
given any v~, a unique / is available for use in the right hand sides 
of the above definitions. 

Using the above definitions, (M) can be expressed using V as 



(1 + a g (t)) — 50'- (/t) = j;(»t) 50 V) +S te (^(,t )>f t )< (72) 



dt 

We now make a technical assumption for the moment, the validity 
of which will be demonstrated later: 

Assumption 2.7: \\£ 50(0 || < M||50(O|| for some < M < °°. 
Assumption 2.7 implies that the magnitude of jjj 50(0 is within a 



constant factor of the magnitude of 50(0, the two are of the same 
order of magnitude. Intuitively, it implies that 50(0 has a bounded 
rate of change. It follows that 50*(r"'") also has a bounded rate of 
change: 

Lemma 2.10: \\ 50"^) || < M2||50 1 "(O)|j for some < M2 < °°. 

6 We will establish later that this assumption is in fact a consequence of the 
external inputs i>*(v) being small. 
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Proof: Using @, @ and Q, we have 



80 1 "(f 1 ^ 



df 



80(f) 



1+Z(f) 

Since |A:(f)| < 1 from Assumption 16, we have S^ft) || < 



m||^ 80(f) || for some m < °°. Using Lemma [110| , the result follows 



(f 1 ") in (ff 



From Lemma p.lOt it is apparent that the term a g -%f 
is of second order, hence can be dropped from the linearization. As 
a result, (^) becomes 



df 



S0 t (f h 



(73) 



Note that (^) is an LPTV system with period T* 



K. Floquet solution of time-shifted LPTV system 
Since |73| ) is > LPTV, the Floquet expressions in Section II-D apply, 
with f and 80(f) replaced by f ' and 80^ (r ), respectively, and f5 ext (f) 
replaced by (0*(ft),ft). (Q becomes 

sp{t) = I MtW^>$ (4) • ^ + (4) + 



scalar 



£ 2,-(' + ) /' ^^vfW-g; (^(T),T)dT, 

1=1 



(74) 



then the solution of the CPS can be expressed as in (^lj), i.e., as 

where <j)*(t) is the periodic, synchronized solution of the externally- 
unperturbed system of oscillators. 50(f), the deviations from the 
orbit of the externally-unperturbed system, remain small and 
bounded for all f (with respect to t he e xternal perturbations {a,(f)}). 

Proof: Subject to Assumption 2.5 and Assumpti on J2.6| , Lemma 
p.8| establishes that bfa(<j)*(fi),t) is small; applying (]70p shows that 
(•,•)> which appears in (Q), is also small. 

Choosing a g (t) to satisfy ( plf ) ensures that cj(f') ( |75| ) vanishes, as 
demonstrated above. As a result, the f = 1 terms in ( |74| ) (which 
correspond to Floquet multiplier p\ = 1 or equivalently, Flo quet 
exponent /ii = 1) remain bounded and small. From Assumption E^l 
and Lemma 2.4, the remaining Floquet exponents /i2, • ■ • all have 
strictly negative real parts. With (•,•) small as noted above, this 
implies that the terms corresponding to i = 2, ■ ■ ■ jiV in (j74j) also 
remain bounded and small for all f. As a result, 80^ (f*) remains 
bounded and small for all f. Applying (jo9[), 80(f) also remains small 
and bounded for all time. This immediately validates Assumption 2.5 
That b$(-, •) is s mall (by assumption) and ( p7| ) holds also validates 
Assumption 2,6l 

Differentiating ([74|) and proceeding in a similar manner, Assumption 
2.7 can also be shown to be valid. 



while the term c\ (f ) in (p^), which causes unbounded growth and 
resulting breakdown of linearization, becomes 



T), t) rfT. 



(75) 



L. Choosing 0C g (t) to circumvent breakdown of linearization 

To avoid unbounded growth of 80^ (f'), which would invalidate the 
present time -shifted linearization procedure in the same manner as 
Section IL-F and Lemma 15 previously, it is imperative that cj(f^) in 
({75]) remain small and bounded (with respect to bX (•,•)). This can 
be achieved by the simple expedient of requiring that the integrand 
in (|75h vanish, i.e., 



v[(t).S^(0*(t),t) = O,Vt. 



(76) 



Substituting f + for T in @, and applying @, (||), ([«]), and 
(Q), we obtain: 

O = v[(t)-S; 2 (0*(t),t) 

^o=v[(f t ).s; 2 (0*(ft),ft) 

^O = vf(ft).[^(0*(?t),f)-z(f)g 1 ( f t)] 
=^ J S:(f)v[(ft).f7 1 (ft)=vf(ft).^(0*(ft) jf ) (77) 
^K(f)=v[(f + ).^(0*(ft),f) 
^a f (f)=v[(f t )-^(0*(ff),f) 



- a f (f) = v[(f + a ? (f))-fy,(0*(f + a g (f)),f) 



From the considerations of Section |ll-l| through (j77j), we are able to 
prove the following Theorem: 

Theorem 1: Given a system of N coupled oscillators modelled in 
the phase domain by the CPS equations ( |ll| ) and mutually injection 
locked, satisfying Assumption [2. l| and Assumption [2.4|. If the external 
perturbations to the system {a;(f)} ( fl~o| ) (equivalently, bq,(-,-) in (fTl])) 
are small, and if cc g (t) is chosen to satisfy (|77|), ;'.e., 

— a f (f) =v[(f + a f (f))-^(0*(f + a g (f)),f), 



III. Conclusion 

( [77] ) and Theorem |l| establish that a g (t), the time shift (or phase) 
of the coupled PPV system, obeys a relationship identical in form 
to the PPV equation (El) for individual oscillators. In other words, 
groups of synchronized oscillators may be abstracted by a single 
"effective PRC/PPV" function that dictates the group's "effective 
phase response" tt g (t) to external perturbations via the single, scalar 
differential equation (|77|). As such, it provides a rigorous basis for the 
empirical practice of measuring PRCs of complex oscillatory systems 
that are synchronized (e.g., [n]). 

( [77] ) may be used to analyze the dynamics (e.g., noise and lock- 
ing/pulling behaviour) of groups of synchronized oscillators, just as 
(g) is used for individual oscillators. Moreover, Theorem [j] may be 
applied repeatedly to abstract the effective PRC/PPV of groups of 
synchronized oscillators that are organized hierarchically over multi- 
ple levels, enabling the development of computationally efficient and 
scalable methods for analysing and abstracting the phase dynamics 
of large systems of coupled oscillators. 
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